Spatio-temporal dynamics of wormlike micelles under shear 
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Velocity profiles in a wormlike micelle solution (CTAB in D2O) are recorded using ultrasound 
every 2 s after a step-like shear rate into the shear-banding regime. The stress relaxation occurs 
over more than six hours and corresponds to the very slow nucleation and growth of the high-shear 
band. Moreover, oscillations of the interface position with a period of about 50 s are observed during 
the growth process. Strong wall slip, metastable states and transient nucleation of three-band flows 
are also reported and discussed in light of previous experiments and theoretical models. 

PACS numbers: 83.60.-a, 83.80.Qr, 47.50.+d, 43.58.+Z 



Wormlike micelle solutions exhibit strong shear- 
thinning behavior due to the coupling between their mi- 
crostructure and the flow. Along the steady-state flow 
curve (shear stress a vs. shear rate 7), a drop of up 
to three orders of magnitude in the effective viscosity 
r/ = a/j is observed in a very narrow stress range lead- 
ing to a stress plateau (for a review, see Refs. 0,0)- 
Such a sharp transition made wormlike micelle solutions 
appear as a model system to study shear-induced effects 
in complex fluids. 

A partial understanding of the flow curve of wormlike 
micelles has emerged thanks to local scattering and ve- 
locimetry experiments. Above a critical shear rate 71, 
a birefringent band normal to the velocity gradient oc- 
cupies an increasing part of the gap as the shear rate is 
increased 0. This shear- induced birefringent band 
was clearly identified as a nematic phase in the case 
of concentrated solutions and has a low viscosity com- 
pared to the isotropic phase 0, 0. Moreover, the ve- 
locity field was shown to separate into two differently 
sheared bands 00: a weakly sheared region flowing at 
71 and a highly sheared region at 72, the upper limit 
of the stress plateau. In the well-documented system 
CPCl/NaSal in brine 0, the equilibrium position of the 
interface between the two shear bands increases linearly 
with 7 consistently with the "lever rule" 0. However, 
such a simple shear-banding scenario remains controver- 
sial due to (i) temporal fluctuations of the flow field and 
(ii) slow transients that question the existence of truly 
stationary states 0, E3 ■ 

At this stage, due to the limited temporal resolution of 
current local velocimetry techniques (Nuclear Magnetic 
Resonance 0,0,0 and Dynamic Light Scattering Q), a 
time-resolved description of the velocity field is still miss- 
ing. In this Letter, velocity profiles obtained in cylindri- 
cal Couette geometry by high-frequency ultrasonic vc- 
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locimetry [ill ] and recorded every 2 s simultaneously to 
global rheological data are presented for a startup ex- 
periment followed over about six hours. For one of the 
much studied micellar systems (CTAB in D2O), we show 
that two very different time scales are involved in the dy- 
namics of shear banding, (i) A three-step nucleation and 
growth of a high-shear band occurs on a few hours, (ii) 
"Fast" oscillations of the band position with a period of 
about 50 s are pointed out. Subtle additional effects such 
as strong wall slip and transient nucleation of three-band 
flows are also reported and discussed. 
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FIG. 1: (a) "Quasistatic" flow curve cr(j) obtained at T = 
44° C under controlled shear rate (15 min per data point) by 
averaging the shear stress over the last 100 s of each step. 
Inset: enlargement over 7 = 60-1000 s _1 showing the stress 
plateau. The dotted line corresponds to the stress response 
to a step-like shear rate from to 200 s~ . o and x symbols 
indicate respectively the maximum and the last value of a. 
(b) Temporal response of the shear stress a(t) after a step- 
like shear rate from to 200 s _1 (thin line). The thick dotted 
line is the sigmoidal relaxation with = 56.4 Pa, <tm = 109 
Pa, and t 3S = 15450 s (see text). Inset: enlargement over 
t — 0-30 s showing the stress overshoot. 



We focus on the salt-free wormlike micelle solution 
made of Cetyl Trimethyl Ammonium Bromide (CTAB) 
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at 20% wt. in deuterated water (D2O) at a temperature 
of 44°C i.e. in the vicinity of the isotropic-nematic (I-N) 
transition that occurs at T IN = 39°C 0,0- Rheological 
data are measured using a stress imposed rheometer (TA 
Instruments AR 1000) and a Couette cell of inner radius 
Ri = 24 mm and gap width e = 1.1 mm. The cell is sur- 
rounded by a solvent trap containing water to prevent 
evaporation. Our local velocimetry technique is based 
on time-domain cross-correlation of high-frequency ul- 
trasonic signals backscattered by the moving fluid. Post- 
processing of acoustic data allows us to record a velocity 
profile in 0.02-2 s with a spatial resolution of 40 /im (see 
Ref. 0] for more details). In order to enhance the scat- 
tering properties of our system, we add a small amount of 
colloidal particles (1% wt. polystyrene spheres of diam- 
eter 3 to 10 ^m). We checked that the linear rheological 
properties as well as the plateau behavior were not sig- 
nificantly affected by the addition of those scatterers. 

Using the shear rate imposed mode of the rheometer, 
we first apply increasing shear rates during a scanning 
time of 900 s per step. The resulting flow curve a(j) 
shown in Fig^a) presents a stress plateau at a* as 84 Pa 
that extends from 71 as 50 s _1 to 72 as 500 s _1 corre- 
sponding to a drop in the effective viscosity by an order 
of magnitude. Note that, in the system under study, the 
features of the plateau (a*, 71, and 72) depend dramati- 
cally on the scanning time as already reported in Ref. Q . 
The bump in (7(7) revealed by the inset of Fig. CJa) is 
characteristic of underlying metastable states |l2j. Thus, 
between 71 and 72 and even for scanning times of 15 min, 
the flow curve of Fig. ^a) does not represent true equilib- 
rium states but rather describes qualitatively the plateau 
behavior of our system. 

In the following, we apply a step-like shear rate into 
the plateau region from to 200 s _1 at time t = and 
we record simultaneously the shear stress response o~(t) 
and the velocity profiles for about six hours. As seen 
in Fig. |Tp3). after an elastic overshoot during the very 
first seconds, a(t) increases from 60 to 105 Pa in 2000 s, 
then remains nearly constant at 107 Pa for 2500 s, and 
finally decreases very slowly for the rest of the exper- 
iment. This ultraslow decrease may be well fitted by a 
sigmoidal function a (t) = <Too + (&m~ &00) exp((— i/r ss ) ) 
with t ss = 15450 s and a ^ = 56.4 Pa [13 ■ Note that 
such a long equilibration time may be surprising for a 
wormlike micelle solution but is consistent with previous 
rheological data on the same salt-free system close to the 
I-N transition We checked that this very long time 
scale cannot be attributed to any evaporation effect by 
measuring the weight fraction of the sample before and 
after shearing using thcrmogravimctric analysis. 

Figure [5] gathers the full ultrasonic velocimetry data 
on a spatio-temporal diagram of the local shear rate 
j(x,t) calculated from the velocity v(x,t) as j(x,t) = 
— (R\ + x) v ^^ x ■ The abscissae correspond to time t 
and the ordinates to the radial position x inside the gap 
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FIG. 2: Local shear rate y(x,t). A linear gray scale is used: 
black and white correspond respectively to 7 = s _1 and 
7 > 300 s _1 . The thick dashed line is given by 5(t) = 5 (to) + 
C /* (ff(t') - <Too)dt', with to = 13500 s, 6 (to) = 0.17 mm, 

? = 1.17 10 -6 mmPa _1 s _1 . a 
from Fig. |Hb). 



of the Couette cell with x — (x = e) at the rotor (sta- 
tor). This representation reveals that the flow field slowly 
evolves in time but also undergoes "fast" fluctuations 
[l3j . By first looking at the velocity profiles averaged 
over at most 6 min, we may distinguish three different 
regimes (indicated by arrows on Fig. |3 and summarized 
in Fig. 0) along the slow stress relaxation. 
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FIG. 3: Description of the slow dynamics. Velocity profiles 
< v(x) > averaged over (a) t = 2-25 (•), 4400-4700 (o), (b) 
12000-12230 (■), and 23600-23840 s (□). The error bars are 
the standard deviations of these estimates. The rotor velocity 
at 7 = 200 s _1 is «o = 196 mms -1 . 

At the inception of the flow and for about 200 s, two 
bands bearing different shear rates, 7^ as 360 s _1 and 
7 as 16 s" 1 , coexist in the gap of the Couette cell (see 
• symbols in Fig.|3Ia)). The interface between the shear 
bands is located at S as 0.3 mm from the stator. More- 
over, wall slip is present since the fluid velocity neither 
reaches the rotor velocity vo = 196 mms -1 at x = nor 
perfectly vanishes at x = e. We define the slip velocities 
at the rotor v s \ and at the stator v S 2 as the differences 
between the fluid velocity near the wall and the corre- 
sponding wall velocity, and get v s i as 75 mms -1 and 
v S 2 as 3 mms . Slip velocities are thus strongly dissym- 
metric and demonstrate that wall slip occurs preferen- 
tially in the low viscosity phase. 

Second, after this first regime, the highly sheared band 
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disappears so that the bulk flow, when averaged over 
a few minutes, becomes homogeneous for about 8000 s. 
This regime is reminiscent of the "milky" turbid phase 
observed at "intermediate" times in Ref. Q . Strong wall 
slip is detected at the rotor (v s i rts 134 mms -1 ) while 
slip at the stator remains negligible (see o symbols in 
Fig- Eta)). Note that the sliding layer at the rotor sup- 
ports a very high shear rate of the order of 3000 s _1 
and thus unambiguously differs from a band of nematic 
phase thinner than our resolution of 40 (im (that would 
flow at 7/j). Another striking feature of this homogeneous 
velocity profile is the value of the "true" shear rate (cal- 
culated by removing the contributions of sliding layers) 
7truc ~ 63 s" 1 , which falls between the local shear rates 
jh and 7; measured previously in the high- and low-shear 
bands. For such a shear rate, the flow is expected to show 
(at least) two shear bands. This suggests that during this 
part of the experiment, the system explores a metastable 
branch of the flow curve located above the stress plateau 
at a* . Indeed, this second stage roughly corresponds 
to the increase and stagnation of the stress response a(t) 
well above both the asymptotic value <7oo = 56.4 Pa and 
the "quasistatic" plateau value er* w 84 Pa (see Fig. 

Finally, from t f=a 8000 s until the end of the experi- 
ment, the velocity profiles show the growth of a low vis- 
cosity layer in the vicinity of the rotor (see Fig. Elb)) 
qualitatively similar to the "long" time behavior reported 
in flow birefringence experiments Q . This highly sheared 
band expands until it fills roughly one third of the gap. 
During this growth stage, the slip velocity v s \ decreases 
from 134 to 64 mms -1 so that the local shear rate in the 
high-shear band 7^ remains roughly constant and equal 
to 360 s _1 . Moreover, the weakly sheared band bears a 
constant local shear rate ji w 30 s _1 . Thus, although the 
asymptotic state is still not perfectly reached, the inter- 
face position at the end of the experiment is compatible 
with the lever rule: S = (jt IUC -7;)/ (% - 7i) e ~ 0.3 mm. 

Let us now get a closer look at the fast dynamics su- 
perimposed to the slow evolution described above. In the 
regime where the flow appears homogeneous on average 
(t w 200-8000 s), Fig. Efa) and (b) reveal that a thin 
band of width 6 < 150 /im is actually nucleated and de- 
stroyed within short intervals of a few seconds unevenly 
every 5 to 20 s. When present, the band is sheared at jh 
and can be clearly distinguished from wall slip. Ref. |4j 
also reports the presence of a small shear band at the 
rotor as the "milky" phase disappears during the second 
regime. 

Moreover, Fig. 0Jc) and (d) show that, during the ex- 
pansion of the highly sheared band (t f=a 8000-24000 s) , 
the interface is in fact subjected to periodic oscillations 
with a period of 50 s. As Fig. ^points out, the slip veloc- 
ity v s i(t) oscillates in phase opposition to the interface 
position 8(i), whereas the true shear rate jtrue(t) is syn- 
chronized with 5(t). These oscillations sometimes give 
way to the nucleation and destruction of a second highly 








i! 






11! 


jj" '*j f||||| i|u 






















I 


1 « mi 





4400 4500 4600 4700 

().4[| i TT l p l WPI : | 



18650 18750 18850 18950 




0.5 

x (mm) 



221(10 



22400 22700 
t(s) 



FIG. 4: Description of the fast dynamics, (a) Velocity profiles 
during the intermittent apparition of a highly sheared band 
at the rotor (x < 150 um) at t = 4505 (•), and 4508 s (o). 
(b) Enlargement of 7(2:, t) over t = 4400-4700 s. (c) Velocity 
profiles during one interface oscillation at t = 18800 (o), 18810 
(•), and 18818 s (□). (d) Enlargement of j(x,t) over t = 
18650-18950 s. (e) Velocity profiles during the nucleation of 
a second highly sheared band at t = 22292 (o), 22350 (•), 
22388 (□), and 22463 s (■). (f) Enlargement of j(x,t) over 
t = 22100-22700 s. 
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FIG. 5: Temporal fluctuations of (a) the position of the inter- 
face S(t), (b) the slip velocities at the rotor v s i(t) (solid line) 
and at the stator v S 2(t) (dotted line), and (c) the true shear 
rate 7true(i), recorded during the oscillations of Fig. Pfd) 
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sheared band at the stator (see Fig.^Jc) and (f)). This 
phenomenon develops within roughly 200 s and takes 
place rather intermittently. According to Ref. 01 1 such 
a transient event is conceivable in Couette geometries 
of low enough curvature. On the spatio-temporal data 
7(x, t), it shows up as a white patch for x > 0.8 mm. 
This process seems to occur more frequently at the end 
of the experiment (see Fig. |5J). Oscillations similar to 
those of Fig. 0Jd) are clearly visible before and after the 
event of Fig.Eff). 

Let us conclude by discussing the main question raised 
by our study: what is the origin of the various time scales 
involved in the spatio-temporal dynamics of our micellar 
system? Concerning the slow time scales, the present 
results reveal a strong similarity with the evolution re- 
ported in Ref. However, our system is more than two 
orders of magnitude slower: "intermediate" and "long" 
times correspond respectively to 200-8000 s and 8000- 
24000 s here vs. 5-35 s and 35-175 s in Ref. 0. Such 
a huge discrepancy could be explained by the differences 
in composition, electrostatic screening, or proximity to 
the I-N transition. In the framework of theoretical ap- 
proaches derived from the Johnson-Segalman model |14j . 
such long time scales mean that the diffusion coefficient 
D of the shear stress through the interface is vanish- 
ingly small. The interface speed is predicted to scale 
as VD(<7 — (Too)- The dashed line in Fig. [21 shows that 
our data are fairly compatible with this prediction. 

Finally the presence of fast time scales ranging from 5 
to 50 s and thus strongly exceeding the micelle relaxation 
time Tfl w 40 ms remains very puzzling |l5l | . Recent phe- 
nomenological models including dynamical equations for 
the micelle length and/or relaxation time predict oscil- 
lations and even chaotic-like states of the flow field that 
could be connected to the dynamics observed experimen- 
tally [lfij . However, we believe that wall slip dynamics 
may play a major role in our experiments. Indeed, time 
scales of the same order have been found after small shear 
rate jumps between two banded states They were 

related to the interface reconstruction after band destabi- 
lization. Here, due to wall slip, the true shear rate is not 
stationary (see Fig. Etc)). The interface could thus be 
periodically destabilized and reconstructed by the same 
process as in Ref. 0|. As a cause for the slip dynam- 
ics, one may invoke some multi-valucdncss of the shear 
rate inside the lubricating layers at a given shear stress, 
which would lead to an instabi lity near the walls similar 
to stick-slip or "spurt" effect [l8j. Further studies will 
focus on trying to isolate the bulk behavior from that of 
sliding layers and determine how the two are coupled. 

In summary, ultrasonic vclocimetry allows the study 
of sheared complex fluids on time scales of about 1 s. 
Such a time resolution is crucial since a wrong picture of 
the flow may be deduced from velocity profiles averaged 



on longer times with slower techniques. In the case of 
a wormlike micelle solution, various dynamical regimes 
were unveiled: the slow nucleation and growth of a shear 
band as well as more complex faster behaviors. These 
results, together with evidence for strong wall slip dy- 
namics, demonstrate the importance of spatio-temporal 
fluctuations in shear-banded flows. More generally, they 
emphasize the need for time-resolved measurements as 
well as dynamical theoretical approaches in the study of 
complex fluid flows. 
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